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Section  I 


Introduction 


Background 


An  important  characteristic  of  any  inventory  control  system  is  the  probability 
distribution  of  forecast  errors.  Such  a  distribution  quantifies  the  percent  of  time  that 
forecast  demand  is  above  or  below  the  forecasted  values,  and  the  magnitude  and 
frequency  of  these  errors.  Such  information  is  critical  in  constructing  cost  effective- 
safety  levels  for  individual  ihventory  items.  At  present,  the  safety  level  calculations 
utilized  in  the  Economic  Order  Quantity  (EOQ)  Buy  Computation  System  (DOS  2)  assume 
that  errors  in  forecasting  the  total  number  of  units  demanded  in  a  procurement 
leadtime  are  normally  distributed.  In  testing  this  assumption,  however,  Demmy  (1979) 
observed  that  aggregate  forecast  errors  were  described  by  a  highly  skewed  distribution. 
One  explanation  for  this  result  is  that  the  demand  process  itself  is  best  described  by  a 
stationary  skewed  distribution.  For  example,  the  negative  binomial  and  several  of  the 
compound  Poisson  processes  have  this  characteristic.  On  the  other  hand,  Hayya  (1980) 
has  suggested  that  the  observed  error  distribution  may  be  the  result  of  combining  the 
biased  forecast  errors  that  result  when  moving  average,  "straight-line"  projections  are 
used  to  forecast  a  mixture  of  items  that  are  in  fact  subject  to  significant  up  and  down 
trends.  Still  another  explanation  is  that  the  demand  process  is  described  by  a  non¬ 
stationary,  skewed  distribution.  There  may  of  course  be  many  other  explanations  for 
the  observed  aggregate  error  distributions. 


.  I«2 

In  this  paper,  we  investigate  the  characteristics  of  forecast  errors  associated  with 
current  D062  forecasting  methods.  Our  goal  is  to  obtain  useful  descriptions  for  the 
distribution  of  demands  for  EOQ  items  that  will  be  convenient  for  use  in  safety  level 
computations.  The  paper  consists  of  several  major  sections.  In  section  I,  we  describe 
the  current  forecasting  methods  utilized  in  the  D062  system  in  more  detail,  and  we 
review  the  current  assumptions  embedded  in  DQ62  safety  level  computations.  In 
Section  II,  we  describe  the  results  of  several  preliminary  investigations  into  the  nature 
of  EOQ  demand  forecast  errors.  This  section  presents  plots  of  ten  year  demand 
histories  for  representative  EOQ  items,  and  also  presents  the  results  of  preliminary 
studies  of  the  consistancy  of  forecast  errors  across  aircraft,  time  periods,  and  demand 
rate  classes.  In  Section  III,  we  present  refined  models  of  the  distribution  of  EOQ 
demand  forecasting  errors  and  we  present  useful  approximation  formulas.  Finally, 
Section  IV  summarizes  the  major  findings  of  this  report. 

Large  amounts  of  data  were  analyzed  during  this  study.  To  reduce  the  volume  of 
this  specific  report,  a  majority  of  the  background  data  is  presented  in  separately  bound 
Appendices. 

Let  us  now  discuss  the  characteristics  of  current  D062  methods  for  forecasting 


EOQ  demands. 


Current  D062  Forecast  Methods 

AFLCR  57-6  contains  a  detailed  description  of  AFLC  EOQ  demand  forecasting 
methods.  In  this  section,  we  briefly  outline  these  procedures.  Basically,  the  current 
D062  forecast  for  recurring  EOQ  demands  is  based  upon  a  eight  quarter  moving 
average,  with  proportional  adjustments  for  projected  changes  in  associated  aircraft 
flying  programs.  For  example,  Figure  1-1  illustrates  an  item  with  an  average  demand  of 
5  units  per  quarter  during  the  past  S  quarters.  Suppose  future  flying  activity  is 
forecasted  to  remain  at  the  same  rate  as  was  observed  during  the  past  S  quarters.  In 
this  case,  D062  would  forecast  5  units  of  demand  per  quarter  for  each  future  quarter,  as 
illustrated  by  dashed  line  in  Figure  1-1.  On  the  other  hand,  if  future  flying  activity  was 
expected  to  be  20  percent  higher  in  the  next  8  quarters  than  that  experienced  during 
the  previous  8  quarters,  the  forecast  quarterly  demand  would  be  (1.20)  x  (5)  =  6  units 
per  quarter. 

Figure  1-1  also  illustrates  the  calculation  of  a  MAD  for  this  hypothetical  item. 
MAD  stands  for  Mean  Absolute  Deviation,  a  measure  of  the  average  difference  between 
the  8  quarter  average  and  each  of  the  8  data  values  that  were  actually  observed.  If  no 
change  in  flying  program  activity  is  expected,  current  forecast  rules  then  use  this 
historical  MAD  as  an  estimate  of  the  MAD  of  demand  in  each  future  quarter.  On  the 
other  hand,  if  future  flying  activity  is  expected  to  be  either  higher  or  lower  than  that 
experienced  in  the  past  two  years,  the  MAD  associated  with  future  quarterly  demand  is 
set  equal  to 


Future  MAD  =  R*83  x  (Old  MAD) 


where  R  denotes  the  ratio  of  projected  flying  activity  during  the  future  two  years  to 
the  flying  activity  observed  during  the  past  two  years.  This  adjustment  is  based  on 
statistical  relationships  observed  by  Robert  Stevens  of  AFLC/XRS  which  indicates  tnat 
the  standard  deviation  of  demand  per  period  is  proportional  to  mean  demand  activity 
raised  to  the  .85  power. 

The  specific  data  values  used  in  computing  the  8  quarter  moving  average  depend 
upon  the  Supply  Management  Grouping  Code  (SMGC)  associated  with  a  given  EOQ  item. 
Figure  1-2  illustrates  the  basic  rules.  For  example,  items  with  SMGC  codes  of  X  and^P^ 
(i.e.  items  with  projected  demand  rates  less  than  $5,000  per  year)  utilize  gross  demands 
in  both  the  average  and  MAD  calculations,  where  gross  demands  are  the  sum  of  sales 
demands,  transfer  demands,  Foreign  Military  Sales  (FMS),  and  non-recurring  demands. 
Serviceable  returns  are  not  included  in  this  calculation.  On  the  other  hand,  items  with 
SMGC  codes  of  P  and  M  (i.e.  items  with  annual  demands  greater  than  $5,000  per  year) 
utilize  net  demand  in  the  average  and  MAD  calculations.  In  this  case,  sales,  transfer, 
and  FMS  demands  are  summed  and  then  past  serviceable  returns  are  subtracted  from 
this  total.  This  gives  net  quarterly  demand  for  each  of  the  past  8  quarters.  Both 
average  and  MAD  values  are  based  on  net  demands  for  SMGC  P  and  M  items. 


Data  Used  in  this  Study 


To  develop  the  distribution  of  forecast  errors  associated  with  the  above  forecasting 
methods  we  utilized  historical  data  from  several  major  sources.  First,  data  describing 
both  forecasted  and  observed  aircraft  flying  programs  were  obtained  from  die  K004  and 
G0333  systems,  respectively,  while  information  on  actual  DOS 2  EOQ  demand  histories 
were  obtained  from  the  INSSIM  Data  Bank.  The  DOS  2  System  maintains  a  weapon  code 
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FORECAST  METHOD  DEPENDS  UPON  SMGC 


SMGC 


FORECAST*  IS  BASED  ON 


X  and  T 


GROSS  DEMANDS 
(sales  *  transfer  ♦  FMS 
♦  non-recurring) 


P  and  M 


NET  DEMANDS 

(sales  ♦  transfer  ♦  FMS 
less  serviceable  returns) 


*Both  Mean  and  MAD  values  ere  computed  using  these  rules. 


Figure  1-2.  Forecast  Calculations  by  Supply  Management 
Grouping  Code. 
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for  each  EOQ  item.  This  code  identifies  the  primary  aircraft  supported  by  the  specific 
EOQ  item.  Actual  aircraft  flying  programs  associated  with  a  given  aircraft  code  for 
the  period  CY73-CY79  were  obtained  through  the  G0333  system  from  output  product  A- 
G033J-PAR-M1-MMO.  This  information  was  collected  for  the  interval  Duly  1970 
through  3une  1979.  Predictions  of  flying  program  activity  which  were  made  for  each 
quarter  throughout  the  1970s  were  obtained  from  the  K004  data  system,  using  report 
K004.D81A  RCS:  NR-LOG-LR(AR)7208.  Reference  2  presents  the  detailed  data  for 
actual  and  predicted  program  activity  collected  in  this  process.  In  addition,  Reference 
2  presents  plots  of  both  the  actual  and  predicted  program  activity,  and  plots  of  the 
ratio  predicted  to  actual  flying  programs  for  each  aircraft  associated  with  the  INSSIM 
Data  Bank.  We  refer  to  this  latter  ratio  as  the  "forecast  accuracy  ratio." 

D062  Demand  Data 


The  item  data  used  in  this  study  were  obtained  from  the  INSSIM  Data  Bank.  This 
data  bank  currently  contains  demand  histories  for  thousands  of  D062  items  managed  by 
the  Sacramento,  Oklahoma  City  and  Warner  Robins  Air  logistics  Centers  (ALC)  for  the 
interval  CY  71-79,  a  total  of  38  quarters  of  demand  data. 

The  INSSIM  Data  Bank  was  constructed  over  a  period  of  years  by  AFLC/XRS  by 
Mr.  Fred  Conway  and  Mr.  Armin  Rubbert.  In  building  this  data  set,  several  problems 
were  encountered  regarding  the  availability  of  data  in  past  years  and  the  readability  of 
the  "old"  tapes  in  the  data  bank.  As  a  result,  the  following  rules  were  adopted  in 
building  the  historical  data  files: 


1.  Sales,  Transfer,  and  FM*>  demands  were  combined  to  provide  a  single  demand 
total  per  quarter. 

2.  Non-recurring  demands  were  not  available  and,  thus,  were  not  included  in  the 
data  bank  records. 

3.  Sales  returns  and  Transfer  returns  were  combined  to  provide  a  single  value  for 
serviceable  returns  Jt>y  quarter. 

4.  Lead  time  and  inventory  management  codes  were  obtained  from  the  records 
corresponding  to  the  first  quarter  of  FY75.  This  includes  the  weapon  code,  the  code 
which  identifies  the  primary  aircraft  associated  with  each  EOQ  item. 

5.  On-hand  and  on-order  assets  were  obtained  from  the  oldest  available  D062 
stock  status  record.  This  corresponded  to  the  first  quarter  of  FY74. 

6.  Items  with  Special  Codes  of  C,  D,  E,  I,  M,  X,  U,  or  N  in  any  fiscal  year  were 
deleted  from  the  file.  Requirements  for  items  with  these  codes  are  computed  using 
manual  methods,  and  consequently  were  not  candidates  for  inclusion  in  the  INSSIM  Data 
Bank. 

7.  Items  with  incomplete  demand  histories  were  also  deleted  from  the  Data  Bank. 
That  is,  an  item  was  included  in  the  INSSIM  Data  Bank  only  if  demand  history  records 
were  present  in  the  D062  system  for  each  of  the  fiscal  years  in  the  interval  CY71-79. 


Hence,  this  rule  eliminates  from  consideration  all  items  which  either  enter  or  leave  a 
specific  ALC's  data  files  during  the  CY71-79  interval,  or  which  were  transferred  from 
one  ALC  to  another  during  the  interval.  Items  which  were  switched  from  management 
under  the  D062  computation  system  to  the  D041  system,  or  vice  versa,  would  also  be 
eliminated  by  this  rule. 


Statistical  S tudies 


Once  the  historical  data  files  were  available,  a  large  number  of  statistical  studies 
were  performed  with  this  data  base.  Figure  1-3  illustrates  the  major  types  of  analyses 
which  were  done.  As  illustrated  in  the  Figure,  our  statistical  analysis  data  system  had 
four  major  components.  First,  as  illustrated  in  Block  A,  we  developed  detailed  item 
plots  of  demands  and  returns  for  several  hundred  specific  items  from  the  INSS1M  data 
bank,  and  we  then,  manually  reviewed  each  of  these  plots  in  an  attempt  to  discover 
underlying  trends  or  other  patterns  that  were  shared  among  the  items.  Reference  3 
presents  plots  of  a  sample  of  one  hundred  items  studied  in  this  process.  As  shown  in 
Block  B,  a  second  major  step  involved  the  analysis  of  standardized  forecasting  errors  Zi 
associated  with  individual  period  forecasts.  In  these  studies,  we  attempted  to  evaluate 
the  consistency  of  forecasting  errors  from  time  period  to  time  period  and  to  discover  if 
there  were  any  significant  error  differences  among  the  distributions  of  period  forecast 
errors  among  aircraft,  time  periods,  demand  classes,  or  other  stratifications  of  D062 
items.  Analysis  of  variance  (ANOVA),  correlation,  frequency  distribution  analysis,  and 
cross  tabulations  were  the  primary  analysis  tools  used  in  this  step.  As  shown  in  Block 
C,  similar  analyses  were  performed  for  the  distribution  of  errors  LTZi  associated  with 
forecasts  for  a  given  number  of  time  periods.  The  results  of  the  analysis  steps 
illustrated  in  Blocks  A,  B,  and  C  are  discussed  in  Section  II.  Based  on  the  results  of 
these  analyses,  we  then  sought  appropriate  transformations  of  the  observed  distribu¬ 
tions  of  forecast  errors  which  could  be  used  to  simplify  the  calculation  of  individual 
safety  level  policies.  This  is  illustrated  by  Block  D  in  Figure  1-3.  The  results  of  these 
efforts  are  discussed  in  Section  III.  Let  us  now  discuss  the  precise  definitions  of  the 
statistics  Zi  and  LTZi  which  were  the  subject  of  these  analyses  efforts. 
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Definition  oi  Standardized  Errors 

In  performing  these  studies  it  was  necessary  to  normalize  our  data  so  that  errors 
associated  with  different  items  might  be  compared.  Two  types  of  standardized  errors 
were  computed.  We  refer  to  these  as  *' period  errors",  which  we  denote  by  the  symbol 
"Zi",  and  the  "leadtime  demand  error",  which  we  denote  by  the  symbol  "LTZi".  The 
computation  formulas  for  these  errors  are  presented  in  Figure  1-4.  As  shown  in  the 
figure,  the  standardized  period  error  Zi  associated  with  quarter  i  is  computed  by 
subtracting  the  forecasted  demand  for  quarter  i  from  the  actual  demand  observed  in 
this  quarter.  This  result  is  then  divided  by  the  8  quarter  historical  MAD.  The 
standardized  leadtime  error  LTZi  is  computed  in  a  similar  manner,  but  in  this  case  total 
forecasted  net  demand  for  the  given  leadtime  is  subtracted  from  the  total  observed 
demand  for  this  interval.  The  result  is  then  divided  by  MAD. 

Numbering  Conventions 


As  noted  above,  38  quarters  of  data  were  available  for  this  study,  covering  the 
interval  CY7 1-CY79.  For  convenience,  we  refer  to  the  first  quarter  of  CY7 1  as  quarter 
1,  the  next  quarter  as  quarter  2,  and  so  on.  This  numbering  scheme  is  illustrated  in 
Figure  3.  Hence,  quarter  9  refers  to  the  first  quarter  CY73  and  quarter  20  refers  to  the 
last  quarter  of  CY75. 
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Figure  1-4.  Formulas  for  Standardized  Errors. 


Since  our  data  base  covered  a  9  year  interval,  a  total  of  4  independent  S  quarter 
base  periods  were  available  in  performing  our  forecast  error  studies.  These  base 
periods  are  the  intervals  CY71-72,  CY73-74,  CY75-76,  and  CY77-78,  respectively.  Two 
of  these  base  periods  are  illustrated  Figure  1-5.  In  performing  our  error  analyses,  we 
began  with  the  CY71-72  base  period.  For  a  given  item,  we  computed  the  average  and 
MAD  of  demands  associated  with  this  base  period.  Next,  we  used  the  appropriate  flying 
program  ratio  to  adjust  this  average  of  CY71-72  demands  to  obtain  a  forecast  of 
quarterly  demands  for  each  quarter  in  the  CY73-CY74  interval.  That  is,  the  data 
observed  in  quarters  1-8  were  used  to  forecast  demands  for  each  quarter  in  the  interval 
from  quarter  9  through  quarter  16.  We  then  used  the  formulas  illustrated  in  Figure  1-4 
to  compute  the  standardized  errors  associated  with  these  forecasts.  The  subscript  i 
associated  with  the  standardized  error  Zi  refers  to  a  forecast  of  i  periods  into  the 
future.  Consequently,  when  the  time  interval  CY71-72  is  the  base  period,  Zl  refers  to 
the  forecast  error  associated  with  a  quarter  9  forecast;  i.  e.  a  forecast  1  time  period 
into  the  future.  Similarly,  Z2  refers  to  the  standardized  period  error  associated  with  a 
forecast  for  quarter  10,  a  forecast  of  2  quarters  into  the  future.  Similar  definitions 
apply  to  the  symbols  LTZi. 

Once  the  standardized  errors  associated  with  each  lead  time  i  were  computed  for 
the  first  base  period,  the  standardized  error  calculations  were  repeated  using  the  period 
CY73-74  as  the  base  period.  In  this  second  case,  the  data  from  quarters  9-16  were  used 
to  forecast  the  quarterly  demands  for  quarters  17-24.  Again  the  standardized  errors  Zi 
were  computed  and  tabulated.  In  this  case,  Zl  refers  to  the  forecast  for  quarter  17, 
while  Z2  refers  to  the  standardized  error  for  the  quarter  18  forecast.  Once  the 
CY73-74  base  period  calculations  were  completed,  similar  calculations  were  performed 


and  tabulated  for  the  base  periods  CY75-76  and  CY77-78.  Detailed  plots  and 
historgrams  for  the  standardized  errors  were  then  produced  from  these  calculations  and 
other  forms  of  statistical  analyses  of  the  standardized  errors  were  performed.  Similar 
activities  were  performed  in  analyzing  the  lead  time  demand  forecast  errors  LTZi. 

The  Current  D062  Model  for  Forecast  Errors 

At  present,  D062  safety  level  calculations  are  based  upon  the  assumption  that  the 
distribution  of  forecast  errors  for  demands  in  a  leadtime  is  normally  distributed.  This 
assumption  is  based  on  the  fact  that  if  a  targe  number  of  identically  distributed  random 
variables  are  summed,  the  sum  will  tend  to  a  normal  distribution  regardless  of  the 
underlying  distribution  of  the  individual  random  variables.  This  is  a  very  common 
assumption  in  inventory  management  systems,  and  simulation  studies  using  actual  D062 
demands  have  shown  that  the  safety  level  formulas  based  upon  these  assumptions 
provide  significantly  more  cost  effective  policies  than  formulas  previously  used  in  the 
D062  system.  Rather  than  use  the  normal  probability  distribution  directly,  however, 
D062  safety  level  calculations  utilize  the  Laplace  distribution  to  approximate  the 
normal.  Presutti  and  Trepp  (1970)  observed  that  the  unit  normal  and  unit  Laplace 
distributions  are  very  close  approximations  to  one  another,  and  they  illustrate  this  fact 
by  comparing  the  cumulative  probability  distributions  shown  in  Table  1-1.  Figure  1-6 
illustrates  the  probability  density  function  for  the  Laplace  distribution.  Notice  that  for 
this  distribution  most  observations  are  concentrated  about  the  mean,  but  the  distribu¬ 
tion  has  a  fairly  long  tail.  Figure  1-7  compares  the  Laplace  distribution  with  the  normal 
probability  distribution.  Note  that  when  compared  to  the  Normal  distribution,  the 
Laplace  has  more  probability  in  the  tails  of  the  distribution  as  well  as  a  slightly  higher 
probability  that  observed  demands  will  be  very  close  to  the  mean. 


Figure  1-7.  A  comparison  of  the  unit  normal  and  unit 
Laplace  (y=0,o«l). 


Figure  1-8  compares  the  cumulative  distribution  function  for  the  normal  probability 
density  with  the  observed  Zi  statistics  from  a  Sacramento  ALC  item  sample.  As  shown 
in  the  figure,  the  cumulative  normal  distribution  will  pass  through  the  fiftieth 
percentile  when  the  standardized  error  Z  =  0.  Also  observe  that  is  extremely  unlikely 
for  a  standardized  error  to  exceed  4  MADS  above  the  mean  of  the  normal  distribution. 
On  the  other  hand,  Figure  1-8  also  illustrates  the  distribution  of  period  errors  Zi 
observed  by  Demmy  (1979).  This  figure  presents  the  observed  distributions  of  Zi  for 
period  errors  in  making  forecast  which  are  1,  2,  3,  and  4  quarters  in  the  future.  These 
particular  curves  were  derived  from  an  analysis  of  a  low  activity  SA-ALC  demand 
sample  utilizing  CY1971  to  72  is  the  base  year  for  the  forecast  calculations.  However, 
almost  identical  curves  were  observed  using  other  Air  Logistics  Center  samples  and 
other  base  periods. 

If  demand  per  quarter  were  in  fact  normally  distributed,  the  observed  distribution 
of  forecast  errors  should  exactly  correspond  with  the  theoretical  normal  curve 
illustrated  in  Figure  1-8.  Obviously,  there  is  a  significant  discrepancy.  One  explanation 
for  this  discrepancy  is  that  the  demands  per  period  is  in  fact  described  by  a  stationary, 
but  skewed  probability  distribution.  Another  explanation  is  that  the  underlying  demand 
process  is  in  fact  stationary  but  subject  to  substantial  trends,  with  some  items  having 
positive  or  increasing  trends  and  other  items  having  negative  (decreasing)  trends.  In 
this  case,  a  "straight-line"  forecasting  technique  would  result  in  a  distribution  of 
forecast  errors  similar  to  that  presented  in  Figure  1-8.  This  explanation  was  suggested 
by  Hayya  (1980).  Still  another  situation  that  would  produce  Zi  curves  similar  to  those 
observed  is  that  demand  is  a  skewed  but  non-stationary  distribution. 


The  major  objective  of  the  studies  reported  in  this  paper  was  to  obtain  a  better 
understanding  of  the  true  nature  of  the  underlying  standardized  error  distributions  for 
Zi  and  LTZi.  In  the  next  section,  we  report  the  results  of  the  preliminary  data  analysis 
required  to  obtain  a  better  understanding  of  the  empirical  data. 


Section  II 


Initial  Results 


At  present,  D062  safety  level  calculations  assume  that  the 
same  probability  distribution-- the  Laplace--is  a  useful  approxi¬ 
mation  to  the  distribution  of  forecast  errors  LTZi  for  all  EOQ 
items.  To  test  the  accuracy  of  this  assumption,  we  performed 
a  number  of  statistical  studies  using  the  D062  demand  histories 
from  the  INSSIM  Data  Bank.  We  were  particularly  interested  in 
determining  how  the  distribution  of  forecast  errors  LTZi  differed 
among  time  periods  and  item  groupings.  For  example,  were  the 
forecast  error  distributions  observed  in  the  early  1970s  similiar 
to  those  observed  during  the  later  part  of  the  decade?  Are  there 
significant  differences  in  forecast  accuracy  across  aircraft 
types?  Or  are  there  significant  differences  in  forecast  accuracy 
across  other  types  of  item  groupings.  If  no  differences  exist, 
a  single  probability  distribution  is  appropriate  for  the  perfor¬ 
mance  of  safety  level  calculations,  and  the  large  amount  of  histo¬ 
rical  data  present  in  the  INSSIM  data  bank  may  be  used  to  develop 
a  precise  description  of  this  distribution.  On  the  other  hand, 
if  significant  differences  exist,  improved  safety  level  calcula¬ 
tions  might  be  developed  by  using  the  specific  probability  distri¬ 
bution  which  is  appropriate  for  the  different  item  classes. 


We  began  our  study  by  developing  plots  of  the  actual  demand 
histories  of  several  hundred  items  from  the  INSSIM  data  bank. 

We  then  visually  compared  these  patterns  both  across  items  and 
with  the  flying  programs  of  the  primary  aircraft  associated  with 
each  item.  In  doing  this,  we  were  primarily  interested  in  obtain¬ 
ing  answers  to  the  following  types  of  questions: 

1.  What  do  D062  item  demand  histories  look  like? 

2.  Are  there  noticeable  trends  in  these  histories?  Are  there 
noticeable  cycles  or  other  patterns  present? 

3.  How  does  this  data  relate  to  flying  hour  programs. 

Decisions  Systems  Working  Paper  81-01  (1981)  presents  plots 
for  100  Sacramento- ALC  items  that  were  studied  in  this  phase, 
while  Figures  I I - 1  thru  11-11  illustrate  our  results.  Figure 
II -1  presents  the  actual  and  predicted  flying  programs  for  F/FB 
111  aircraft  for  the  CY71-79  interval,  while  Figures  II-2  thru 
II-ll  present  representative  demand  histories  of  items  which 
support  the  F/FB  111.  Figures  1 1  -  2  thru  1 1  -  5  are  '’high”  demand 
items;  that  is,  these  items  had  annual  dollar  demand  rates  which 
exceeded  $5,000  per  year  in  the  CY  71-72  interval,  while 
Figures  II-6  thru  II-ll  present  plots  for  items  which  had  demands 
of  less  than  $5,000  per  year  during  this  interval. 


Figure  II-2.  Observed  Demands  and  Returns  for  FSN1560  000528153. 


u'ini'iN  %  <v .  <  >> 


i  ft  «  •  <  *  •  »  •  #., 


• t>  uc 

s:x  'r 

-•  1>  »w 

«•  u>  -j 


CO  «v  •  wl  *v  •  •»  • 


« 

Ol  O'  M 

¥- 

•  «l  * 

« 

ft.  ✓  Al 

O 

4DAAUJMV  m 

«. 

Ml 

v> 

»- 

« 

r» 

r» 

4 

ft  4  M 
•  m>M  n%  nu> 

fk’ClSlfCS 

y 

u 

AJ 


i-miAAlMi 

•U  •*-« 


.  Figure  II-3.  Historical  Demands  and  Returns  for  a  High  Activity  F/FB111  Item 
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Figure  II-5.  Historical  Demands  and  Returns  for  a  High  Activity  F/FB111  Item 


Figure  II-7.  Historical  Demands  and  Returns  for  a  Low  Activity  F/FBlll  Item 
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As  shown  in  II-l,  the  F/FB111  was  being  introduced  into  the 
Air  Force  inventory  during  the  beginning  of  the  1970s.  The  fly¬ 
ing  program  for  this  aircraft  increased  during  the  first  two 
and  a  half  years  of  the  decade ,  and  has  been  undergoing  a  slight 
downtrend  since  that  time.  As  shown  on  the  right  hand  side  of 
Figure  II-l,  Air  Force  G033J  forecasts  of  flying  activity  for 
the  F/FB111  have  tended  to  exceed  the  actual  hours  flown. 

Figure  1 1 - 2  illustrates  the  observed  demands  and  returns 
for  FSN  1560  000528153,  an  item  which  supports  the  F/FB  111  air¬ 
craft.  In  this  figure,  the  solid  line  denotes  demands  observed 
each  quarter  during  the  CY71-79  interval,  while  the  dotted  line 
denotes  serviceable  returns  recorded  by  the  D062  system.  For 
this  item,  demands  appear  fairly  stable  throughout  the  nine 
year  time  period.  However,  no  noticeable  relationship  with  the 
actual  flying  program  for  the  F/FB  111  is  present. 

Figure  II-3  illustrates  another  type  of  pattern  observed 
in  D062  demand  histories.  If  one  takes  a  large  time  period  view 
one  might  conclude'  that  this  pattern  represents  a  stationary 
but  highly  erratic  demand  pattern.  On  the  other  hand,  if  one 
assumes  that  the  underlying  demand  process  is  fairly  stable  but 
has  a  trend  that  shifts  through  time,  one  might  conclude  that 
this  patten  is  associated  with  an  item  which  initially  has  a 
significant  increasing  trend  during  the  first  three  years  of 
the  time  interval,  followed  by  a  significant  decreasing  trend 
during  the  next  three  year  interval.  The  final  demand  values 
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might  be  interpreted  as  fairly  stationary  demands  with  a  spike 
in  quarter  28,  --  i.e.,  a  large  unexpected  demand  inconsistent 
with  the  underlying  pattern.  For  this  item,  observe  that  if 
one  assumes  trends  are  present  in  the  data,  these  trends  bear 
no  clear  relationship  with  the  flying  program  activities  for 
the  F/FB111 . 

Figure  1 1 - 4  illustrates  another  pattern  found  in  D062  demand 
histories.  This  item  experiences  a  significant  downtrend  during 
the  CY71-75  interval.  It  is  followed  by  nine  quarters  of  zero 
demand  and  then  another  time  interval  in  which  there  are  demands 
in  three  of  the  six  subsequent  quarters.  Again,  the  demand 
histories  for  this  item  bear  no  clear  relationship  to  the  flying 
program  for  the  F/FB111, 

Figures  I I - 5  presents  still  another  pattern.  In  this  case 
demands  are  quite  erratic  with  a  long  term  increasing  trend; 
however,  a  sharp  drop  in  the  level  of  demand  occurs  in  the  last 
few  quarters  of  the  38  quarter  interval. 

Figures  I I - 6  thru  11-11  present  representative  demand  his¬ 
tories  for  ’’low”  activity  items,  i,e.,  items  which  had  demands 
of  $5,000  per  year  or  less  during  the  CY71-72  interval.  The 
vast  majority  of  D062  items  fall  into  this  catagory.  Observe 
that  demands  for  these  items  are  much  more  erratic- -and  therefore 
much  more  difficult  to  predict-- than  the  demands  for  the  high 
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activity  items  discussed  above.  In  general,  demands  for  these 
items  contain  a  large  number  of  periods  with  no  observed  demands. 
When  demands  do  occur,  they  often  involve  quantities  signifi¬ 
cantly  different  than  one  unit.  For  example.  Figure  1 1 - 7  shows 
only  two  quarters  out  of  thirty-eight  in  which  demands  are 
greater  than  zero.  One  of  these  demands  was  for  60  units.  We 
hypothesize  this  pattern  was  due  to  a  requisition  which  was 
improperly  coded,  for  notice  that  60  units  were  returned  one  year 
after  the  60  units  of  demand  were  placed  on  the  D062  system. 

After  reviewing  a  very  large  number  of  plots  for  D062  items, 
it  was  difficult  to  distinguish  any  single  pattern  which  appeared 
useful  in  forecasting.  For  some  items,  demand  does  in  fact 
appear  to  follow  flying  program  activity.  For  other  items, 
apparent  trends  in  demand  appear  to  be  unrelated  to  program  activ 
ity.  For  still  other  items,  no  consistent  trends  were  observed. 
The  only  pattern  which  appeared  common  across  large  numbers  of 
D062  items  was  the  presence  of  spikes-- i.e. ,  the  presence  of 
very  large  and  unexpected  demands  which  are  inconsistent  with 
demand  activity  both  before  and  after  the  occurance  of  the  spike. 
Spikes  are  particularly  bothersome  for  low  activity  items.  Unfor 
tunately,  spikes  appear  in  the  patterns  of  a  large  number  of 
D062  demand  histories. 


I 


The  Distribution  of  Period  Errors  Zi 


In  this  section,  we  discuss  our  preliminary  analysis  of  the 
distribution  of  period  errors  Zi.  Of  particularly  interest  in 
this  section  are  the  following  questions: 

Are  there  significant  differences  in  the  distribution  of 
period  errors  Zi 

(a)  across  base  periods? 

(b)  across  lead  times? 

(c)  across  Air  Logistics  Centers? 

(d)  across  dollar  demand  classess? 

To  answer  these  questions,  we  built  four  item  samples  which 
we  denote  as  OC.H,  OC.L,  SM.H,  and  SM.L,  respectively.  Sample 
OC.H  represents  a  sample  of  INSSIM  Data  Bank  high  activity  items 
managed  by  Oklahoma  City  Air  Logistics  Center,  while  sample  OC.L 
represents  a  sample  of  low  activity  Oklahoma  City  D06Z  items. 
Similarily,  samples  SM.H  and  SM.L  were  obtained  from  the  INSSIM 
Data  Bank  histories  for  Sacremento  Air  Logistics  Center.  To 
be  included  in  a  "high"  activity  category,  an  item  must  have 
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had  an  average  of  at  least  $5,000  of  demands  per  year  during 
the  CY71-72  inteval;  otherwise,  the  item  was  classed  in  the  low 
demand  category. 

As  noted  in  section  I,  a  total  of  38  quarters  of  data  were 
available  for  the  study.  This  permitted  us  to  compute  the  forecast 
errors  associated  with  four  independent  two-year  base  periods. 

The  base  periods  involved  were  the  intervals  CY71-72,  CY73-74, 
CY75-76,  and  CY77-78.  As  discussed  in  chapter  I,  we  used  the 
demands  observed  in  a  given  base  period  to  compute  the  mean  and 
MAD  associated  with  these  demands.  This  information  was  modified 
for  forecasted  changes  in  flying  activity  to  obtain  the  forecast 
demand  in  each  future  quarter.  The  standardized  error  Zi  represents 
the  ratio  (observed  demand  -  forecasted  demand)/MAD,  where  the 
observed  demand  is  associated  with  the  period  i  quarters  beyond 
the  last  quarter  of  the  base  period.  For  example,  in  obtaining 
forecast  associated  with  the  base  period  CY71-72,  the  symbol 
Z2  represents  the  standarized  error  associated  with  demands  observed 
in  the  second  period  of  CY73. 

To  obtain  an  understanding  of  the  shape  of  the  distribution 
of  forecast  errors  Zi,  we  computed  the  standardized  errors  Zi 
associated  with  each  item  sample  and  each  demand  period,  and 
for  forecasts  of  i*l,2,...,12  quarters  in  the  future.  Plots 
of  the  cumulative  distributions  associated  with  each  of  the  error 
statistics  Zi  were  then  obtained  and  studied  in  detail.  In  addi¬ 
tion,  we  used  the  error  statistics  Zi  to  perform  Analysis  of 


Variance  (ANOVA)  studies  to  identify  differences  among  different 
item  groupings.  These  ANOVA  studies  sought  to  determine  if  there 
are  any  significant  differences  in  the  individual  Zi  statistics 
across  aircraft,  base  periods,  or  demand  classes. 

The  detailed  ANOVA  tables  obtained  from  these  studies  are 
presented  in  Appendix  A. 

We  found  significant  differences  in  the  error  characteristics 
associated  the  CY75-76  base  period  relative  to  the  other  base 
intervals.  In  addition,  error  statistics  associated  with  F104 
and  F5  aircraft  appeared  to  be  of  significantly  greater  magnitude 
than  those  associated  with  the  other  twenty-one  aircraft  in  the 
INSSIM  data  bank.  A  possible  explanation  for  these  significant 
differences  is  perhaps  explained  by  the  data  presented  in  Figures 
11-12  and  11-13.  Figure  11-12  presents  actual  and  forecasted 
flying  for  the  F104  during  CY71-79  interval.  Note  that  this 
weapon  system  phased  out  of  the  Air  Force  inventory  about  the 
middle  of  the  period,  but  that  foreign  military  sales  resulted 
in  increased  programs  during  the  latter  part  of  the  1970's. 

Notice  that  particularly  severe  errors  are  associated  with  the 
forecast  program  accuracy  in  quarter  twenty.  This  particularly 
bad  forecast  of  future  flying  program  activity  would  significantly 
bias  forecasts  associated  with  that  period.  Now  observe  the 
predicted  and  actual  flying  programs  for  the  F5  presented  Figure 
11-13.  As  shown  in  the  figure,  the  F5  experienced  a  continuously 
increasing  program  throughout  the  1970's.  The  FS  is  the  only 
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aircraft  in  the  INSSIM  data  bank  which  has  an  ascending  program. 

(Note:  the  F/FB111  has  an  increasing  program  during  CY71-72, 
but  declines  gently  during  the  remainder  of  the  70s).  Errors 
in  forecasting  future  flying  program  activity  for  the  F5  are 
also  very  significant  during  the  75-76  time  period. 

Based  on  the  above  data,  it  appeared  that  demand  data  associ¬ 
ated  with  F-104  and  F-5  aircraft  would  differ  significantly  from 
that  associated  with  other  aircraft.  Consequently,  all  items 
associated  with  these  aircraft  were  deleted  from  future  forecast 
studies.  We  then  repeated  the  Analysis  of  Variances  studies 
for  differences  across,  aircraft  and  across  base  periods,  but 
this  time  we  excluded  F-104  and  F-5  items.  With  these  items 
removed,  there  were  very  small  differences  among  the  mean  forecast 
errors  both  by  aircraft  and  base  period.  Since  thousands  of 
forecast  were  involved,  these  small  differences  were  still  statis¬ 
tically  significant,  but  they  were  of  such  small  magnitude  as 
to  have  little  practical  usefullness.  For  example,  the  proportion 
of  variation  in  mean  forecast  errors  explained  by  differences 
across  aircraft  and  base  periods  are  generally  much  less  than 
one  percent  for  all  ^i  values  studied  and  never  exceeds  two-percent. 
Hence,  it  appears  that  differentiating  with  respect  to  either 
aircraft  or  base  periods  would  provide  little  improvement  in 
estimating  mean  forecast  errors. 

We  also  performed  Analysis  of  Variance  studies  to  determine 
if  there  were  any  significant  differences  in  forecast  errors 


based  upon  the  level  of  activity  observed  during  the  base  period. 

The  detailed  ANOVA  tables  produced  in  this  effort  are  presented 
in  Appendix  A.  We  found,  yes,  there  are  significant  differences 
in  forecast  errors  based  upon  the  observed  base  period  demands. 

If  little  or  no  demand  is  observed  in  the  base  period,  actual 
demands  which  occur  in  future  periods  tend  to  be  greater  than 
the  forecast.  On  the  other  hand,  if  very  high  levels  of  demand 
are  observed  in  the  base  period,  future  demands  tend  to  be  less 
than  the  forecast.  In  fact,  it  appeared  that  the  higher  the 
demand  forecast  (i.e.  the  higher  the  demand  observed  in  the  based 
period),  the  greater  the  tendency  to  overforecast  future  demands. 

In  section  III,  we  will  futher  investigate  differences  among 
forecast  error  distributions  based  upon  demand  classes. 

The  Analysis  of  Variance  studies  indicate  that  little  improve¬ 
ment  in  the  accuracy  of  forecasting  mean  forecast  errors  would 
be  obtained  by  differentiating  across  aircraft  or  base  periods. 

None -the -less,  the  shapes  of  the  error  distributions  might  still 
be  significantly  different  across  these  item  groupings.  To  investi¬ 
gate  this  issue,  we  again  developed  the  empirical  distributions 
of  period  errors  Zi  associated  with  each  item  sample  and  each 
base  period.  This  time,  however,  all  items  associated  with  F- 
104  and  F5  aircraft  were  deleted.  Tables  II-l  and  I I - 4  and  Figures 
11-14  through  11-21  illustrate  our  results  for  the  high  activity 
sample  from  Oklahoma  City  ALC.  Similar  data  for  the  OC.L,  SM.H, 
and  SM.L  samples  are  presented  in  Appendix  B. 
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Table  II-l.  Cumulative  ProbabM  i  t ies  for  Ep^ogram  Factors 

for  1971-72  Base  Year  Forecasts  With  General  Progra 
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Table  1 1 - 2 .  Cumulative  Probabilities  for  Period  Errors  Z 

for  1973-74  Base  Year  Forecasts  With  General  Program  Factor 

for  Sample  OC.H 
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Table  IX- 3.  Cumulative  Probabilities  for  *  Factor! 

for  1975-76  Base  Year  Forecasts  With  General  Program 

for  Sample  OC.H 
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Table  II-4.  Cumulative  Probabilities  for  Period  Errors  Z 

for  1877-78  Base  Year  Forecasts  With  General  Program  Factors 

for  Sample  OC.H 
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Let  us  now  consider  Table  I I - 1 -  Column  1  of  this  table  presents 
the  Cumulative  Distribution  Function  (CDF)  associated  with  fore¬ 
cast  errors  Zl,  i.e.,  the  observed  standardized  error  Zi  when 
data  for  base  year  CY71-72  is  used  to  forecast  demands  for  the 
first  quarter  of  CY73.  Similarly,  Column  2  of  this  table  presents 
the  CDF  of  forecast  errors  associated  with  Z2,  i.e.  errors  in 
forecasting  demands  during  the  second  quarter  of  CY73.  Finally, 
column  12  denotes  the  CDF  associated  with  Z12,  the  distribution 
of  forecast  errors  when  CY71-72  data  is  used  to  forecast  demands 
twelve  quarters  into  the  future  .  Note  that  Samples  OC.H  contained 
a  total  of  3153  items.  The  mean,  variance,  coefficient  of  varia¬ 
tion,  skewness,  and  kurtosis  statistics  associated  with  each 
of  these  distributions  is  presented  at  the  bottom  of  the  table. 
Observe  that  the  variance  of  each  of  these  distributions  is  extremely 
high,  that  the  distributions  tend  to  be  skewed,  and  that  the 
kurtosis  values  are  particularly  large.  In  contrast,  if  the 
CDF  of  forecast  errors  were  normally  distributed,  we  would  expect 
to  see  mean,  standard  deviation,  skewness,  and  kurtosis  values 
of  0,  1,  0,  and  3,  respectively. 

Tables  1 1 - 2  through  1 1 - 4  present  similar  information  for 
the  CY73-74 ,  CY75-76,  and  CY77-78  base  year  forecast,  respectively. 

Figures  11-14  thru  11-21  present  graphs  of  the  information 
presented  in  Tables  I I - 1  through  II-4.  Figure  11-14  presents 
plots  of  period  errors  Zl,  Z2,  Z3,  and  Z4  associated  with  CY71- 
72  base  year  forecasts  for  sample  OC.H.  As  shown  in  the  figure. 
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these  CDFs  appear  to  be  identical.  Similarly,  Figure  11-15  plots 
the  cumulative  distribution  function  for  the  period  errors  Z4, 

Z6,  Z8 ,  Z10,  and  Z12  for  CY71-72  forecasts  for  sample  OC.H. 

Note  that  the  CDFs  for  these  Zi  appear  to  be  very  similar  to 
the  Zl  through  Z4  curves.  Figures  11-16  through  11-21  present 
similar  plots  associated  with  the  other  base  year  forecasts. 
Observe  that  there  is  very  little  difference  among  these  curves 
across  all  base  periods  and  lead  times. 

As  noted  above,  similar  information  is  presented  in  Appendix 
B  for  the  OC.L,  SM.H,  and  SM.L  samples. 

Gene ral  Observation s 


After  analyzing  this  information  in  detail,  we  believe  that 
several  general  conclusions  may  be  reached.  First,  within  a 
given  item  sample  and  base  period,  the  distribution  of  period 
errors  Zi  appear  to  be  almost  identical  for  all  i.  There  is 
a  minor  tendency  to  over-forecast  more  as  the  forecast  horizon 
increases,  but  this  tendency  appears  slight.  Second,  for  all 
items  and  samples,  all  of  the  curves  are  very  similar.  This 
is  particularly  true  for  Zi  values  that  are  greater  than  zero. 
Finally, the  low  demand  item  samples  appear  to  have  a  slightly 
higher  probability  of  over-forecasting  than  for  the  higher  activity 
samples.  In  the  next  section,  we  perform  studies  to  identify 
the  impact  of  differing  demand  rates  upon  the  standardized  error 
distributions. 


Section  III 


Analysis  of  Forecast  Errors  for  Demand  in  a  Leadtime 
for  High  Activity  Items 


In  Section  ny  we  observed  significant  differences  among  item  groups  based  on  average 
dollar  demands  in  a  base  period.  This  is  not  surprising.  Anaylsis  of  item  demand  histories 
shows  that  many  low  demand  D062  items  have  many  periods  of  zero  demands,  but  that 
when  demands  do  occur,  they  are  often  of  large  magnitude  relative  to  the  average  usage 
rate.  On  the  other  hand,  high  activity  D062  items  appear  to  be  less  erratic  though  still 
highly  variable.  As  a  result,  we  conducted  two  basic  categories  of  studies:  (1)  analyses  of 
items  with  relatively  high  usage  rates  of  three  units  per  quarter  or  more,  and  (2)  studies 
of  items  with  lower  levels  of  activity.  In  this  Section,  we  consider  the  high  usage  rate 
items,  while  Section  IV  presents  our  results  for  low  activity  items. 


Definition  of  Standarized  Forecast  Errors 


We  are  interested  in  identifying  the  characteristics  of  errors  in  forecasting  the  total 
units  demanded  in  a  given  number  of  time  periods.  From  probability  theory,  it  is  well 
known  that  if  demands  per  period  are  identically  distributed  random  variables,  then  the 
expected  value  of  a  sum  of  t  of  these  variables  equals  t  times  the  expected  value  of 
demand  in  a  single  period.  Similarly,  if  the  variables  are  independent,  the  variance  of  the 
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sum  is  equal  to  t  times  the  variance  of  demand  in  an  individual  period.  In  symbols,  let  Xj 
denote  the  demand  in  period  i,  and  suppose  Xj  has  mean  u  and  standard  deviation  sx  for  all 
periods  i.  Let  D  denote  the  total  demand  observed  in  t  periods;  i.e. 

(1)  D  =  X,  ♦  Xj  ♦  . .  .  ♦  X, 

Since  the  X.  are  independent  and  identically  distributed, 

(2)  E  (D)  =  t  •  E  (X.)  =  t  u 
and 

(3)  Var  (D)  =  fVariX.)  =  t  sx2 

Hence,  the  standard  deviation  of  O  is  given  by 

sd  ■ 

From  the  studies  discussed  in  Section  II,  we  know  that  period  errors  are  correlated. 
Consequently  we  do  not  necessarily  expect  relationship  (3)  above  to  hold.  Nevertheless, 
we  found  it  useful  to  standardize  forecast  errors  in  a  manner  similar  to  that  employed 
when  demands  are  uncorrelated.  Specifically,  let  R  denote  the  forecasted  quarterly 
demand  rate  associated  with  each  future  period,  and  let  MAD  denote  the  Mean  Absolute 
Deviation  associated  with  the  observed  demand  in  the  base  period.  Then  we  define  the 
standardized  forecast  errors  for  demand  in  the  lead  time  as: 

[Total  Demand  Observed  j 

in  t  quarters  I  ^ * 


MAD  v'T 
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or. 


CZt  =  D  -  R  t 

mad  yr 

If  forecast  errors  are  in  fact  normally  distributed  —  as  assumed  in  the  current  D062 
saftey  level  calculations  —  then  the  observed  distribution  for  the  standardized  error  CZ 
should  be  normally  distributed  with  a  mean  of  zero  and  a  standard  deviation  of  1.25  (since 
for  the  standard  normal  distribution  the  MAD  is  approximately  equal  to  .8s). 

Demand  Class  Numbering  Comentions 

As  noted  above,  we  observed  significant  differences  in  the  forecast  error  characteris¬ 
tics  among  differing  demand  classes.  Consequently,  we  assigned  individual  item  forecasts 
to  specific  demand  classes.  Table  111-1  presents  the  assignment  of  Histogram  Numbers  to 
specific  forecast  demand  rate  categories.  For  example,  histograms  1  thru  6  were 
associated  with  items  which  had  demand  forecasts  in  the  range  of  3  to  6  units  per  quarter. 
Histogram  1  developed  the  cumulative  distribution  function  (CDF)  associated  with  CZj, 
i.e.  the  distribution  of  errors  in  forecasting  demands  for  one  quarter  into  the  future. 
Histogram  two,  on  the  other  hand,  tabulated  ttte  histogram  of  standardized  errors  CZ2 
associated  with  errors  in  forecasting  the  total  demand  in  quarters  1  and  2  combined. 
Similarly,  tables  3,  4,  5  and  6  tabulated  the  CDFs  associated  with  standardized  errors  in 
forecasting  the  total  demands  in  3,  4,  5,  and  6  quarters,  respectively.  In  all  of  these 
cases,  the  forecast  demand  rate  was  in  the  3  to  6  units  range.  Similarity,  Histograms  7 
thru  12  were  assigned  to  items  with  forecast  demand  rates  in  the  range  of  6  to  10  units 
per  quarter.  For  items  in  this  demand  rate  category,  table  7  tabulated  the  CDF  of 
standardized  errors  CZj  for  forecasts  one  quarter  into  the  future,  Histogram  8  tabulated 


Table  III-l 


Assignment  of  Histogram  Numbers  in  Analysis 
of  CZj/k  Distributions 


Histogram 


Forecast 


Numbers 


Demand  per  Quarter  (Units) 


1-6 

7-12 

13-18 

19-24 

25-30 

31-36 

37-42 


3-6 

6-10 

10-31 

31-100 

100-310 

310-1000 

1000-3100 
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the  errors  in  forecasting  total  demands  for  two  future  quarters,  and  so  on.  Other 
histogram  number  assignments  and  the  corresponding  forecast  demand  rate  categories  are 
illustrated  in  Table  HI-1. 

Let  k  =  MAD  /T7  We  calculated  the  statistic  CZ/k  for  each  of  the  four  item  samples 
SM.L,  SM.H,  OC.L,  and  OC.H,  and  tabulated  our  results  using  the  Histogram  numbering 
scheme  shown  in  Table  HI-1.  Our  results  for  Sample  SM.L  are  presented  in  Table  III-2. 
For  example,  consider  the  second  and  third  columns  of  the  table.  Column  2  presents 
values  for  the  CZ/k  statistic,  while  column  3  presents  percentage  counts  for  histogram  1, 
i.e.  counts  for  forecasts  in  the  3-6  units/qtr  range.  As  shown  in  these  columns,  CZ/k 
values  were  less  than  or  equal  to  -1  in  26%  of  the  forecasts  in  this  range,  76%  of  the  CZ/k 
values  were  less  than  or  equal  to  0,  and  86%  of  the  CZ/k  values  were  less  than  or  equal  to 
1.  Other  percentage  values  for  CZ/k  in  the  range  from  -9  to  10  may  be  read  from  column 
3.  Other  columns  of  Table  III-2  present  similar  statistics  for  each  of  the  Histograms 
defined  in  Table  01-1. 

Since  a  picture  is  often  worth  one  thousand  words,  let  us  now  consider  the  shapes  of  these 
distributions. 


Observed  CZ/k  CDFs  for  Sample  SM.L 
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Plots  of  CZ/K 


Let  us  now  consider  Figure  I I I - 1 .  The  curves  shown  on  the 
left  hand  side  of  this  figure  denote  the  cumulative  Distribution 
Functions  (CDF)  of  CZ/K  associated  with  items  having  a  forecast 
demand  in  the  range  of  three  to  six  units  per  quarter,  while 
the  plots  on  the  right  hand  side  of  this  figure  denotes  analy¬ 
tical  approximations  to  this  curve.  As  shown  in  the  figure, 
for  values  of  N  greater  than  zero,  the  observed  CDFs  are  almost 
identical  for  all  six  lead  times;  however,  significant  differ¬ 
ences  among  the  curves  exist  for  negative  N  values.  Figure  III- 
1B  presents  three  approximations  to  the  CDFs  presented  in  Figure 
III-1A.  First,  the  solid  line  denotes  the  average  for  all  six 
of  the  curves  presented  in  Figure  III-1A.  The  large  dashed  line 
presented  in  Figure  III-1B  presents  the  cummulative  distribution 
function  associated  with  the  normal  probability  distribution. 

For  the  normal  distribution  the  fiftieth  percentile  occurs  when 
the  random  variable  N  equals  zero,  and  the  99.99  percentile  occurs 
for  N  equals  to  3.62  standard  deviations  or  (3.62)  (1.25)«  4.52 
MADS.  Finally,  the  third  curve  presented  in  Figure  II -IB  is 
extremely  difficult  to  see,  because  it  lies  almost  on  top  of 
the  average  (solid  line)  curve.  This  third  curve- -represented 
by  a  line  of  small  dashes--is  the  curve  of  a  an  exponential  func¬ 
tion  which  passes  through  the  average  curve  at  the  point  where 
N  equals  zero  and  N  equals  five.  Note  that  this  exponential 
approximation  is  an  excellent  fit  to  the  observed  data. 


Figure  III-1.  CZ/k  Curves  for  SM.L  Forecasts  of  3-6  Units  per  Quarter 
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Figures  1 1 1  - 2  thru  1 1 1 -6  presents  results  similar  to  Figure  III- 
1  for  demand  classes  in  which  the  forecast  is  in  the  range  of 
6-10,  10-31,  31-100,  100-310,  and  310-1000  units  per  quarter, 
respectively.  Each  of  these  curves  present  individual  CDFs  associ¬ 
ated  with  the  CZ/K  statistic  for  each  of  the  six  different  lead 
times.  In  addition,  average,  normal,  and  exponential  approxi¬ 
mations  to  each  of  these  curves  are  also  presented.  Note  that 
all  of  these  curves  are  very  similiar.  In  no  case  does  the  average 
CDF  appear  to  be  normally  distributed.  On  the  other  hand,  the 
exponential  approximation  provides  an  excellent  fit  to  the  upper 
tail  of  the  observed  cumulative  distribution  function  in  all  of 
these  figures. 

If  each  of  the  curves  presented  in  Figures  I I I - 1  thru  Hi- 
fi  are  compared  across  demand  rate  groups,  we  observed  that  all 
of  the  curves  are  very  similiar.  For  example,  Table  m-3  presents 
a  comparison  of  the  "average1'  curves  associated  with  each  demand 
rate  group,  and  displays  the  maximum  differences  among  these 
curves.  As  shown  in  the  table,  there  is  practically  no  differ¬ 
ence  between  the  curves  for  N  values  greater  than  zero.  However, 
differences  as  large  as  71  exist  for  N  values  less  than  zero. 
Fortunately,  in  computing  safety  levels,  our  primary  interest 
is  for  N  values  which  are  zero  or  greater.  Consequently,  for 
safety  level  calculation  purposes,  it  is  useful  for  us  to  compute 
an  average  upper  tail  curve  for  this  class  items. 


Figure  III-2.  CZ/k  Curves  for  SM.L  Forecasts  of  6-10  Units  per  Quarter 


Figure  III-3.  CZ/k  Curves  for  SM.L  Forecasts  of  10-31  Units  per  Quarter. 


Curves  for  SM.L  Forecasts  of  31-100  Units  per  Quarter 


TABLE  1 1 1 - 3 

COMPARISON  OF  AVERAGE  CURVES  OF  P[  cz/k  *  z  1 

FOR  SAMPLE  SMIL 


P  £  cz/k  ^  z  ]  by  Units/Qtr  Category 


3-6  6-10  10-31 


MAXIMUM 

31-100  DIFFERENCE 


l  2  3  4  4  2 

i  6 

7 

9 

11 

5 

l  18 

20 

24 

25 

7 

50 

50 

54 

54 

4 
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Comparisons  with  Other  ALC  Samples 

In  the  above  paragraphs,  we  presented  data  associated  with 
the  SM.L  sample.  However,  similiar  findings  apply  to  the  other 
three  samples  used  in  this  study.  The  detailed  plots  associated 
with  each  of  these  other  samples  are  presented  in  Appendix  C. 

In  all  cases,  we  observe  that  an  exponential  function  is  an  excel¬ 
lent  approximation  to  the  average  curves  for  N  values  greater 
than  zero.  Similiarly,  we  also  find  that  a  single  curve  provides 
a  good  approximation  to  the  exponential  tail  across  all  demand 
rate  groups.  Unfortunately,  it  appears  that  each  item  sample 
is  charaterized  by  its  own  exponential  tail.  To  see  this,  let 
us  consider  Table  III-4  and  Figure  III-8.  Table  III-4  presents 
the  exponential  tail  coefficient  estimates  associated  with  all 
samples,  while  Figure  III-8  plots  each  of  these  coefficient  esti¬ 
mates.  As  shown  in  the  table,  the  coefficient  estimates  are 
very  similiar  within  an  item  grouping,  but  differ  significantly 
from  one  item  group  to  another. 

Additional  Characteristics  of  Forecast  Errors 


Figures  III-9  thru  1 1 1 - 1 1  present  several  other  interesting  chara- 
teristics  of  the  distribution  of  forecast  errors.  Figure  III- 
9  plots  the  mean  value  of  CZ/K  by  demand  category  for  the  OC.H 
sample.  Let  us  consider  the  six  points  on  the  left  hand  side 
of  this  figure.  These  points  represent  the  Mean  forecast  of 
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Table  III-4 


Exponential  Coefficient  Estimates  by  Demand  Class  and  Group 


Model:  F  (z) 

«  1 

-  A  exp 

(Bz) 

Forecast 

CZLTOC.L 

CZLTOC.H 

Demand /QTR 

Histogram 

N 

A 

B 

N 

A 

B 

3-6 

1-6 

3401 

.302 

-.381 

1184 

.412  - 

.409 

6-10 

7-12 

2070 

.292 

-.415 

1172 

.407  - 

.433 

10-31 

13-18 

2806 

.307 

-.476 

3000 

.398  - 

.496 

31-100 

19-24 

1204 

.313 

-.481 

2619 

.387  - 

.535 

16-310 

25-30 

300 

.300 

-.599 

1274 

.392  - 

.483 

31-1000 

31-36 

75 

.310 

-.629 

565 

.397  - 

.467 

1000-3100 

37-42 

14 

.262 

-.402 

237 

.422  - 

.471 

3100-10000 

43-48 

2 

- 

- 

120 

.317  - 

.536 

WEIGHTED  AVERAGE 

9876 

.303 

-.436 

10171 

.397  - 

.485 

Forecast 

CZLTSM.L 

CZLTSM.H 

Demand /OTR 

Histogram 

N 

A 

B 

N 

A 

B 

3-6 

1-6 

2013 

253 

-.463 

360 

.313  - 

.469 

6-10 

7-12 

1045 

.245 

-.444 

244 

.345  - 

.439 

10-31 

12-18 

1235 

.255 

-.452 

505 

.310  - 

.566 

31 -100 

19-24 

394 

.253 

-.416 

203 

.327  - 

.640 

100-310 

25-30 

92 

.245 

-.473 

59 

.273  - 

.743 

310-1000 

31-36 

22 

.302 

-.512 

21 

.372  - 

.299 

36-42 

2 

- 

- 

4 

- 

- 

1000-3100 

37-42 

2 

- 

- 

4 

- 

- 

WEIGHTED  AVERAGE 

4801 

.252 

-.453 

1442 

.320  - 

.530 

GRAND  WEIGHTED  AVERAGE 

N«26 , 290 

A*. 331 

B»  - 

.463 

Mea 


CZ/K  associated  with  forecasts  which  are  at  the  three  to  six 
units  per  quarter  forecast  range.  The  six  points  represent  lead 
times  of  one,  two,...,  six  quarters,  respectively.  Note  that 
for  these  points,  the  mean  CZ/K  value  is  positive--that  is,  actual 
demands  tend  to  be  greater  than  the  forecasted  value.  Further, 
this  mean  error  increases  as  the  lead  time  increases.  The  second 
group  of  six  points  corresponds  to  forecasts  in  which  the  fore¬ 
cast  demand  rate  is  in  the  range  of  six  to  ten  units  per  quarter. 
Note  that  a  similiar  phenomena  occurs  here;  however,  the  mean 
CZ/K  values  are  less.  On  the  other  hand,  observe  that  for  the 
high  demand  rate  catagories,  the  mean  CZ/K  values  are  negative 
--indicating  actual  demands  are  on  the  average  less  than  the 
values  forecast.  Also  note  that  this  over-forecasting  charate- 
ristic  becomes  more  severe  as  the  forecast  horizon  lengthens. 

Figure  III-10  plots  the  standard  deviation  of  the  normalized 
value  CZ/K  verses  the  Histogram  number.  Notice  that  the  standard 
deviation  increases  as  the  lead  time  increases,  but  that  the 
standard  deviation  is  similiar  across  all  demand  rate  catagories. 

If  demand  were  in  fact  independently  distributed,  the  standard 
deviation  of  CZ/K  should  not  change  as  a  function  of  lead  time. 
Finally,  Figure  III-ll  presents  a  plot  of  the  squared  value  of 
skewness  verses  the  kurtosis  associated  with  each  of  the  forcast 
catagories.  As  shown  in  the  figure,  items  with  demand  rate  fore¬ 
casts  in  the  3-6  unit  range  have  skewness/kurtosis  character¬ 
istics  which  lie  between  those  associated  with  gamma  and  log-normal 
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Skewness2  Vs  Kurtosis  for 


functions.  This  suggests  that  perhaps  a  gamma  or  log-normal 
distribution  may  be  a  useful  approximation  to  the  actual  CZ/K 
distribution  for  items  in  this  range.  Note  that  as  the  demand 
rate  catagory  increases,  the  resulting  distributions  tend  to 
be  less  skewed  and  the  kurtosis  measures  also  decrease  slightly. 
However,  the  skewness /kurtosis  values  still  depart  significantly 
from  the  values  that  would  be  expected  from  a  normal  distribution 
of  forecast  errors. 


Section  iV 


Analysis  of  the  Distribution  of  Demand  per  Quarter  for 
Low  Activity  Items 

As  noted  in  earlier  sections,  we  observed  significant  differences  in  forecast  error 
distributions  among  items  grouped  on  the  basis  of  demand  activity  in  a  base  period.  Anal¬ 
ysis  of  item  demand  histories  shows  that  many  low  demand  D062  items  have  many  periods 
of  zero  demands  but,  when  demands  do  occur,  they  are  often  of  large  magnitude  relative 
to  the  average  usage  rate.  On  the  other  hand,  high  activity  D062  items  appear  to  be  less 
erratic  though  still  highly  variable.  Consequently,  it  appears  that  an  optimum  D062 
inventory  management  system  might  utilize  different  methods  in  modeling  high  versus  low 
activity  items.  For  example,  in  many  commercial  inventory  systems,  the  Poission 
distribution  is  used  to  model  low  activity  item  demands,  while  the  Normal  distribution  is 
used  to  model  high  activity  demand  processes.  In  Section  III,  we  developed  an  empirical 
model  of  the  probability  distribution  of  forecast  errors  for  high  activity  items.  In  this 
section,  we  report  the  results  of  similar  efforts  to  develop  an  appropriate  model  for  the 
distribution  of  forecast  errors  for  low  activity  items. 

Initial  Studies 

We  began  our  studies  in  this  area  by  asking  several  questions.  In  particular,  we  asked 
"Are  there  significant  differences  in  the  distribution  in  actual  demands  in  a  given  quarter: 

a.  As  the  lead  time  increases? 

b.  As  the  number  of  demands  observed  in  the  last  eight  quarters  increases? 

c.  As  the  forecast  interval  changes  ?" 

To  answer  these  questions,  we  first  developed  histograms  of  the  net  demands  observed  in 
a  given  quarter  given  that  a  total  of  T  demands  have  been  observed  in  the  past  eight  quarters. 
For  example,  suppose  that  a  total  of  T  =  10  demands  were  observed  during  the  past  eight 
quarters  for  a  given  item.  We  then  tabulated  that  number  of  demands  which  occurred  one 
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quarter  in  the  future  for  this  item.  A  similar  tabulation  was  done  for  all  items  in  which  a 
total  of  T  =  10  demands  had  been  observed  in  the  last  eight  quarters.  This  provided  a 
histogram  of  actual  net  demand  one  quarter  in  the  future  given  a  total  of  10  demands  had 
been  observed  in  the  past  eight  quarters.  This  histogram  construction  process  was  then 
repeated  for  forecasts  of  demand  two  quarters  in  the  future,  three  quarters  in  the  future, 
etc.  Let  X.  denote  the  number  of  demands  actually  observed  in  period  i.  Then  the  above 
tabulation  process  provides  an  empirical  measurement  of  P(X./T),  the  conditional  prob¬ 
ability  of  observing  X.  units  in  quarter  i  given  a  total  of  T  units  were  demanded  in  the 
past  eight  quarters. 

To  obtain  answers  to  the  above  questions,  we  first  developed  histograms  of  P(X./T) 
for  all  items  with  total  CY71-72  demands  T  of  1,  2,  4,  and  8  units,  respectively.  For  each 
group  of  items,  we  developed  separate  histograms  for  forecast  leads  times  of  i  -  1,  4,  and 
eight  quarters,  respectively.  This  resulted  in  a  total  of  (4  base  demands)  x  (3  lead  times)  =12 
separate  histograms  for  a  given  item  sample.  This  histogram  building  process  was  then 
repeated  using  the  intervals  CY73-74,  CY75-76,  and  CY77-78  as  the  base  period.  This 
resulted  in  a  total  of  48  separate  P(Xj/T)  histograms  for  each  item  sample. 

We  next  carefully  analyzed  and  compared  these  histograms.  We  found  that  within  a 
given  item  sample,  demand  per  quarter  tends  to  be  identically  distributed  for  a  given  base 
demand  T.  That  is,  the  observed  cumulative  distribution  functions  for  PfX^T)  are  very 
similar  for  all  lead  time  values  i  (i=  1,4,8)  and  for  all  the  base  period  intervals.  However, 
items  with  very  few  base  period  demands  T  have  significantly  different  PCX^/T)  curves 
than  items  with  high  levels  of  past  demand  activity. 

Based  on  the  above  results,  we  developed  another  set  of  histograms.  In  this  effort, 
we  assurmned  that  P(X./T)  was  identically  distributed  across  lead  times  i  and  across 
forecast  intervals.  Histogram  1  was  used  to  tabulate  actual  quarterly  net  demand  for 
items  with  total  (past)  base  period  demands  of  T=1  unit,  while  histogram  two  tabulated  net 
quarterly  demands  for  items  with  past  demands  Ts2  units.  A  total  of  48  histograms, 
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corresponding  to  total  base  period  demands  of  T=i,  2,  ...,  48  units,  respectively,  was 
constructed  ior  each  item  sample. 

We  next  sought  analytical  approximations  to  these  empirical  results. 

The  Poisson  Model 

The  Poisson  probability  distribution  is  often  used  to  model  the  distribution  of 
demands  for  low  activity  items.  Let  R  denote  the  average  quarterly  demand  rate  for  a 
given  item.  Then  if  demands  are  Poisson  distributed,  the  probability  of  x  units  being 
demanded  in  a  given  quarter  is  given  by 

p(x)  =  RX  exp(~R)  x=0,l,2,... 
x! 

If  T  units  have  been  demanded  in  the  past  eight  quarters,  an  estimate  of  R  is  T/S. 

Figures  IV-1  and  IV-2  plot  the  probability  of  zero  demands  in  a  given  quarter  for 
samples  SM.L  and  OC.L  for  both  the  Poisson  distribution  and  for  the  empirical  P(Xj|T) 
histograms.  As  shown  in  the  figures,  observed  D062  demands  have  a  significantly  higher 
probability  of  zero  demands  than  is  predicted  by  the  Poisson  model.  For  the  Poisson 
model,  as  the  demand  rate  increases,  P(XsO)  becomes  very  small.  However,  for  the 
observed  D062  demands,  P(X=0)  was  30%  or  higher  even  when  40-50  units  had  been 
demanded  in  the  past  two  years. 

As  a  result  of  these  observations,  we  concluded  that  the  Poisson  distribution  is  not  a 
good  model  for  describing  the  demand  process  for  low  activity  D062,  items. 

Exponential  and  Split  Poisson  Approximations 

Since  the  observed  data  was  clearly  not  Poisson  distributed,  we  sought  convenient 
mathematical  approximations  to  the  observed  data.  Two  approximation  methods  were 
tested.  First,  we  obtained  an  exponential  approximation  by  fitting  an  exponential  curve 
through  the  P(X=Q)  and  P(X=S)  points  on  the  observed  cumulative  distribution  function  for 
net  demand.  As  a  second  approximation,  we  fit  a  probability  model  that  consisted  of  a 
spike  of  probability  at  XsO  (corresponding  to  zero  demands),  with  the  remaining 
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probability  allocated  according  to  the  Poisson  formula.  We  refer  to  this  second 
approximation  as  the  "split  Poisson,"  since  we  have  split  the  probability  values  between 
the  P(X=0)  point  and  the  Poisson  model. 


Testing  the  Approximation 


To  test  the  above  approximations,  we  developed  a  series  of 
graphs  comparing  the  observed  cumulative  distribution  functions 
of  net  demand  and  the  two  approximations  described  above.  Because 
of  the  limited  amount  of  data  available,  the  CDF  for  any  particular 
PC  X  I  T  )  curve  was  often  rather  unstable.  To  avoid  this  problem, 
we  constructed  an  average  curve  by  pooling  the  data  for  several 
histograms.  We  then  constructed  approximations  to  the  average 
curve. 

Our  results  for  histograms  1  thru  6  are  presented  in  Figure 
IV- 3.  To  the  left  of  this  figure,  we  present  the  observed  cumula¬ 
tive  distribution  functions  for  the  histogram  in  which  the  demands 
observed  during  the  last  eight  quarters  equal  to  one  unit,  two 
units,  ...»  six  units,  respectively.  We  then  computed  the 
average  of  these  six  curves,  and  plotted  the  average  curve  as 
well  as  the  exponential  and  spilt  Poisson  approximations  to  this 
curve.  These  latter  results  are  shown  on  the  right  hand  side 
of  Figure  IV- 3.  As  shown  in  the  figure,  the  exponential  curve 
provides  a  reasonable  approximation  to  the  average  curve,  while 
the  spilt  Poisson  provides  a  less  desirable  fit.  The  coefficients 
associated  with  the  estimation  equations  are  shown  at  the  bottom 
of  this  figure.  The  coefficients  a  and  b  correspond  to  the  expo¬ 
nential  approximation  model  P(X  *  x)  ■  1  -  a  expC  b  x),  while 
the  parameter  c  is  a  parameter  of  the  spilt  Poisson  model. 


Figures  IV-4  thru  IV-10  present  similar  results  for  histograms 
7  thru  48.  In  each  case,  on  the  left  hand  side  we  display  the 
observed  cumulative  distribution  functions  for  net  demand  given 
that  demand  in  the  past  eight  quarters  equals  to  a  specified 
value  T,  while  on  the  right  hand  side  we  compare  the  average 
of  these  six  curves  with  the  exponential  and  spilt  Poisson  approxi¬ 
mations.  Notice  that  is  the  number  of  demands  in  the  base  period 
increases,  the  exponential  approximation  becomes  a  better  fit 
to  the  actual  CDF.  On  the  other  hand,  the  fit  of  the  split  Poisson 
becomes  very  poor  as  the  base  period  demands  increase.  Consequently, 
we  conclude  that  the  exponential  approximation  would  be  a  reason¬ 
able  model  for  demand  in  a  given  quarter  for  low  demand  D062 
items  and  that  the  spilt  Poisson  is  not  useful  for  this  purpose. 

Table  IV-1  presents  the  exponential  coefficient  estimates 
associated  with  the  Oklahoma  City  and  Sacramento  ALC  low  activity 
samples  while  Figure  IV- 11  presents  a  graph  of  this  data.  As 
shown  in  the  figure,  the  coefficient  estimates  increase  as  the 
number  of  units  in  the  past  eight  quarters  increase.  Least  squares- 
fit  equations  which  might  be  used  to  obtain  A  and  B  coefficients 
for  a  given  value  of  base  demand  T  are  presented  at  the  bottom 
of  Table  IV-1. 
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Observed  Cumulative  Distribution  of  Net  Demand 


Observed  Cumulative  Distribution  of  Net  Demand 
for  Sample  SM.l  and  Base  Demand  of  25-30  Unit 
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Figure  IV-10.  Observed  Cumulative  Distribution  of  Net  Demand 
for  Sample  SH.l  and  Base  Demand  of  43-48  Units 
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Exponential  Coefficient  Estimates 

for 

Low  Activity  Items 


Model:  P(Net  Demand  4  x)  ■  1.  -  A*exp(B*x) 


Total  Units 

Demanded  in  _ 

Base  Period  Midpoint  A 


OC.L 
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3.  S 

.19 

-.281 

.187 

7-12 

9.5 

.293 

-.239 

.320' 

13-18 

15.5 

.353 

-.212 

.392 

19-24 

21.5 

.393 

-.184 

.442 

25-30 

27.5 

.420 

-.164 

.460 

31-36 

33.5 

.477 

-.164 

.493 

37-42 

39.5 

.435 

-.158 

.575 

42-48 

45.5 

.497 

-.151 

.600 

Least  Squares  Estimate: 

A*. 2239  ♦  0.006S  T  A-  .2155  +  0.0089  T 

B"-. 266+0. 0029  T  B-  -.3328  +0.0045  T 

where  T*  Midpoint  of  Base  Demand 


Exponential  Fit  Coefficient  Estimates 


SECTION  V 


Summary  of  Results 

This  paper  has  presented  statistical  data  describing  the  distri¬ 
bution  of  forecast  errors  associated  with  current  D062  demand 
forecasting  procedures.  Section  I  provides  a  detailed  review 
of  these  procedures  while  Section  II  describes  the  results  of 
preliminary  studies  of  the  nature  of  EOQ  demand  forecast  errors. 
Finally,  Sections  III  and  IV  present  refined  models  for  the  distri 
bution  of  forecast  errors  for  high  and  low  activity  items,  respec¬ 
tively. 

In  Section  III,  we  studied  forecast  errors  for  items  with  fore- 
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cast  demand  rates  of  three  units  per  quarter  or  more,  while  Sec¬ 
tion  V  describes  the  error  distributions  for  lower  activity 
items.  In  Section  III,  we  found  that,  for  positive  standardized 
forecast  errors  (i.e.,  errors  in  which  demand  exceeds  the  fore¬ 
cast  value),  an  exponential  model  provides  an  excellent  fit  to 
observed  error  statistics.  An  exponential  approximation  may 
also  be  useful  in  approximating  the  distribution  for  negative 
forecast  errors,  but  the  fit  is  not  as  good  for  this  side  of 
the  cumulative  error  distribution  curve.  Note  that  if  leadtimes 
are  considered  fixed  and  safety  levels  are  restricted  to  be  non¬ 
negative,  then  the  positive  error  portion  of  the  error  distri¬ 
bution  function  is  all  that  is  needed  to  set  cost-effective 
safety  levels.  For  this  situation,  we  have  an  excellent  fit. 
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In  Section  IV,  we  also  found  an  useful  exponential  approximation. 
We  defined  P(X^IT)  as  the  conditional  distribution  of  net  demand 
in  a  given  quarter  i  given  that  a  total  of  T  units  of  demand 
were  observed  in  the  p;ast  eight  quarters.  We  found  that  E(XilT) 
appears  to  be  identically  distributed  across  leadtime  periods 
i  and  across  forecast  intervals.  Further,  we  found  the  model 
P(X*x|T)  *  1  -  a  exp  (bx)  for  X  -  0,  1,  2,  to  be  a  good 
approximation  to  the  observed  P(X^|T)  distributions.  Analytical 
or  numeric  calculation  of  convolutions  of  this  model  might  then 
be  used  to  estimate  the  distribution  of  cumulative  net  demand 
for  any  given  leadtime. 
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Differences  Between  Aircraft  Groups 


Table  A-2.2.  Breakdown  by  Aircraft  Code  for  All  Items  in  Sample  OC.L. 

Analysis  of  Variance  for  Differences  Between  Aircraft  Groups 
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Table  A-3.2.  Breakdown  by  Aircraft  for  Sample  OC.L.  Excluding  F104  §  F5  Items. 
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Table  A-5.3.  Breakdown  by  Dollar  Demand  Class  for  All  Items  in  Sample  SM.H. 
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Table  A- 5. 4.  Breakdown  by  Dollar  Demand  Class  for  All  Items  in  Sample  SM.L. 

Analysis  of  Variance  for  Differences  Between  Dollar  Demand  Classes 
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Empirical  Distribution  of  Period  Errors  Zi 
Excluding  F104  and  F5  Items 
for  Samples  OC.L,  SM.H,  and  SM.L 
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Cumulative  Probabilities  for  Period  Errors  Z 
for  1971*72  Base  Year  Forecasts  With  General  Program  Factor^ 
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